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ABSTRACT 

Aims. We present spectroscopic observations obtained with the infrared Spitzer Space Telescope, which provide insight into the H2 
physics and gas energetics in photodissociation Regions (PDRs) of low to moderate far-ultraviolet (FUV) fields and densities. 
Methods. We analyze data on six well known Galactic PDRs (L1721, California, N7023E, Horsehead, rho Oph, N2023N), sampling 
a poorly explored range of excitation conditions (x ~ 5 - 10^), relevant to the bulk of molecular clouds in galaxies. Spitzer observa- 
tions of Ht rotational lines are complemented with Ht data, including ro-vibrational line measurements, obtained with ground-based 
telescopes and ISO, to constrain the relative contributions of ultraviolet pumping and collisions to the H2 excitation. The data analysis 
is supported by model calculations with the Meudon PDR code. 

Results. The observed column densities of rotationally excited H2 are observed to be much higher than PDR model predictions. In 
the lowest excitation PDRs, the discrepancy between the model and the data is about one order of magnitude for rotational levels J >3. 
We discuss whether an enhancement in the H2 formation rate or a local increase in photoelectric heating, as proposed for brighter 
PDRs in former ISO studies, may improve the data-model comparison. We find that an enhancement in the H2 formation rates reduces 
the discrepancy, but the models still fall short of the data. 

Conclusions. This large disagreement suggests that our understanding of the formation and excitation of H2 and/or of PDRs energet- 
ics is still incomplete. We discuss several explanations, which could be further tested using the Herschel Space Telescope. 
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1. Introduction for this discrepancy (iHabart et al.L |2004 . With H2 UV ab- 

sorption line observations performed by the Far Ultraviolet 
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^ ; The bulk of mterstellar matter is found m regions of low-to- Spectroscopic Explorer (FUSE), the H2 formation and excita 

(N| moderate opacity to UV and visible light where stellar radia- ^^^^ processes in diffuse and translucent clouds have also been 

m ' tion governs the chemical and thermal state of the gas. These reconsidered. In particular, UV data show significant amounts 

in ; photodissociation or photon-domin^d regions (PDRs, for a re- rotationally excited H2, for which UV photons could not 

(N view see lHoUenbach & Tielen4 | 1999|) are responsible forrepro- (jjg u^jq^g Seating source ( Grv et all 120021 iNehme et al.L 

^ cessing much of the energy output from stars, reemitting this ^he H2 excitation in the diffuse interstellar medium 

O energy in the infrared-millimeter wavelengths, including a rich tracing the dissipation of in terstellar turbulence in C- 

^ , mixture of gas lines (e.g., emission in fine-structure, rotational, shocks ('Flower & Pineau D es ForetsL [ToOSh or interstellar vor- 

IJ • and rovibrational lines). They are privileged objects for studying ^j^g^ (Falgarone et al 2005) 

■ the chemical and physical processes of the interstellar- medium There'has been very litde exploration of the physics of PDRs 

(ISM). The motivation of this study is to test our understanding ^-^^ moderate FUV fields and densities, an intermediate regime 

H of the excitation of H2 in a regime of space parameters that has ^^^^^^^ ^j^^ ^-^^^^ ^j^^^^ ^j^^ ^^-^^^^ pi^^^ ^j^^ 

. >^ been poorly studied. ^j^^ j^frared spectrograph (IRS) onboard the Spitzer Space 

Through its ability to observe pure rotational lines of H2 to- Telescope provides a unique opportunity to observe H2 pure 

wards a number of bright PDRs, the Infrared Space Observatory rotational line emission in low-brightness sources. As part of 

(ISO) has given important information on the local micro- jhe SPECPDRQ program dedicated to the study of PDRs with 

physics of interstellar- gas. One conclusion of these studies is Spitzer, we used this capability to study the H2 Une emission in 

that the H2 hne intensities and the gas temperature, derived regions with x ~ 5-10^ times the Solar Neighborhood Far- 

from the first rotational levels of H?,, are higher than rnodel cal- uV interstellar- radiation field as given by .Drain e (1978). Our 

culations (Bertoldi, 1997; Draine & Bertoldi, 1999"; Thi et al., ^^mple consists of well-known PDRs, so that they are good test 

[1999^ Kemper et al., , 1999 ; , Li et al., , 2002; , Habart et al.. , 2003) . f.^^^^ for models. These observations are analyzed in combina- 

Changes in the description of heating and cooling processes 

or of the H2 formation rate have been proposed to account 1 See http://www.cesr.fr/~joblin/SPECPDR public/SPECPDR.html 

' For a discussion on the definitions of the mean inters tellar radia- 

Send offprint requests to: E. Habart. e-mail: emilie.habart@ias.u- tion field used in the liter ature on PDR s see A ppendix B of lAllen et al.l 

psud.fr O004) or Appendix C of iLe Petit et all ( l2006h . 
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tion with previous ISO and ground-based data of rotational and 
rovibrational H2 lines when available, in order to help constrain 
the relative roles of ultraviolet pumping and collisions in estab- 
lishing the level populations. 

The paper is organized as follows. In Sect.|2] we present our 
PDR sample. In Sect. [3] we present the IRS spectrometer obser- 
vations. In Sect, m we briefly describe the PDR model used to 
analyze the data, and in Sect. |5] we compare the observed H2 line 
emission to the PDR model predictions. In Sect.|6] we compare 
our results with UV absorption measurements in diffuse clouds. 
In Sect. |7] we discuss several possibilities of explaining our re- 
sults. In Sect. [8] in light of our PDR observations, we discuss 
H2 line infrared emission measurements in galaxies and along 
Galactic lines of sight. Our conclusions are summarized in Sect. 

El 



2. Selected PDRs 

In the framework of the SPECPDR program, we have observed 
several photodissociation regions - a total of eleven including 
nine galactic and two extra-galactic sources - with all the instru- 
ments onboard Spitzer. In this paper, we select six of the nine 
galactic sources. They are close to the Sun {d ~ 100-600 pc) 
with well-defined geometry and radiation field. Most of them 
are inclined with respect to the line of sight with nearly edge- 
on (shell-like) geometry. They represent classic examples of the 
interaction between UV radiation and an interstellar molecular 
cloud. They span a wide range of physical conditions; tempera- 
ture of the exciting star: T,ff = 15, 000 - 37, 000 K, FUV field: 
X - 5 - 1000 in Draine units, proton gas densities: « = 10^-10^ 
cm""*. A brief description of each source is given below. 

L1721. At the low excitation en d, our sample contains the 
Lynds Dark Nebula 1721 (L1721. ' Lvndi Il962h . a nearby iso- 
lated molecular cloud in the p Ophiuchi region. This moder- 
ate opacity cloud is heated on one side by the star v Sco and 
by a more isotropic interstellar r adiation field create d by the 
Upper Scorpius association (e.g., iHabart et al.L 1200 lb . The B2 
star V Sco is located at ~130 pc from the Sun. Visual extinction 
map and emissions of the dust and gas cooling lines ([C^] 158 
j um, [O^l 63 um) h ave been measured toward the L1721 cloud 
dHabart et al.Ll2001i) . The exciting radiation field is estimated to 
hex =3 to 6 times the Draine radiation field. 



California. The California Nebula (NGC 1499) is an emission 
nebula located in the constellation Perseus. The 07.5III star ^ 
Per is the source of radiation of the California nebula. This star 
is located at ~540 pc. Here, we focus on the filament situated 
at 45' from the star The intensity of the exciting radiation field 
after dilution and assuming that the distance of the exciting star 
to the PDR is equal to the distance projected onto the sky is 
estimated to be;^^ -19 times the Draine field. 

N7023E. NGC 7023 is a reflection nebula around the Herbig 
B5Ve illuminating star HD 200775. This star is located at -430 
pc. Three PDRs are visible around the star Here, we focus on the 
NGC 7023 East located further east of the stai- (170"). Spitzer 
spectroscopic dust studies of the East PDR was performed by 
dBerne et al l l2007h . The exciting radiation field is estimated to 
be 87 times the Habing field, i.e. ;t' ~51 times the Draine field. 

Horsehead. The Horsehead nebula (B33) emerges from the edge 
of the L1630 molecular complex. The cr Ori star, a 09.5V 



system, is the main heating source in this region. This star 
is located at ~350 pc. We focus here on the western edge of 
the nebula situated at 0.5° from cr Ori. Emissions of the gas 
and dust associated with the western edge have been mapped 
from the visible t o millimeter wavelength s by various authors 
(e.g..lPound et all 12003: Tevssie r et ail 120 04: Petv et al., 2004; 
Habart et al.', '2005; Ward-Thompson et al., 2006; Philip p et al . 
2006; Goicoechea et al., 2006; Pety et al., 2007). This' PDR is 
illuminated and viewed nearly edge-on. The intensity of the ex- 
citing radiation field after dilution and assuming that the distance 
of the exciting star to the PDR is equal to the distance projected 
onto the sky is x ~60 times the Draine field. 

Rho Oph. The Ophiuchi molecular cloud is a nearby star- 
forming region. The western part of this cloud is bounded by a 
bright filament structure, forming an edge-on (shell-like) PDR, 
illuminated by the B2 HD 147889 star ( Abergel et al.Lll999 ). The 
HD 147889 star is located at ~ 130 pc. We focus here on the west- 
ern edge situated at ~600" from the illuminated star ISO spec- 
troscopic dust and gas (H2) studies of this PDR were performed 
(Habart et al. , 200 3). The exciting radiation field is estimated to 
hex ~250 times the Draine field. 



N2023N. NGC 2023 is a reflection nebula excited by a B 1.5V 
star (HD37903), which is embedded in the LI 630 cloud. The 
HD37903 star is located at ~470 pc. In this work, we focus on 
the filament situated at 150" north of the star Spitzer spectro- 
scopic dust studies of the N2023 north PDR was performed by 
Compiegne et al. (2008). The intensity of the exciting radiation 
field after dilution and assuming that the distance of the exciting 
star to the PDR is equal to the distance projected onto the sky is 
estimate d to be y ~550 times the Draine field. According to the 
study of Compieg ne et al] ( l2008h . there is a visual extinction of 
Ay equal to 1.25 between the star and the filament. The intensity 
of the exciting radiation field after dilution and extinction is 120 
times the Habing field, i.e.^ ~70 times the Draine field. 

In Table [T] we summarize the distances, radiation sources, 
gas densities and geometries for each region as determined 
from the h terature. The d i stance s have been estimated by 
Hipparcos ('Perrvman et al.L Il997h and other measurements 
(Hoogerwerf et al., 200Q). As given above for most of our 
sources, the radiation field x factor is determined from the ex- 
pected FUV luminosity of the exciting star and assuming that 
the distance of the exciting star to the PDR is equal to the dis- 
tance projected onto the sky. This is an upper limit in principle. 
Uncertainties in x due to the errors in the distance of the excit- 
ing star from the Sun are given, so, estimations of x based on 
observations of dust and gas lines emissions are also reported. 
The gas density for each region is also given in Table [T] The 
density inferred from various atomic/molecular species shows 
a relatively large dispersion. This dispersion could result from 
density gradients from the o uter hot layer to the inner molecular 
cold layer (see, for example. Waknslev et alll2000tlHabart et all 
2003, 2005). However, what is needed for our study is the den- 
sity in the H2-emitting region. 

3. Observations 

3.1. IRS spectrometer data 

The selected PDRs were observed using the IRS spectrometer in 
the low (A/AA -60-127) and/or high (A/AA -600) spectral res- 
olution mode. In the low spectral resolution mode, spectral line 
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Table 1. Sample of PDRs with distance, effective temperature of the exciting star, incident FUV radiation field strength, gas density, 
and geometry. 



Source 


Dist." 


Exciting star 


FUV field* 


Density'^ 


Geometry'' 


Ref. 




(pc) 


(Teff, spec.type) 


(Draine) 


(cm-3) 






L1721 


134!- 


20,000 K, B2V 


4.5 [3-6] 


3 10'-3 10' 


S-C 


(1) 


California 


543+""' 


37,0000, 07.5III 


19 [7-37] 


lO^-lO* 


E-O 


(2) 


N7023E 


15,000 K, B5Ve 


51 [27-82] 


7 10^ 


E-O 


(3) 


Horsehead 


352+1^* 


33,000 K, 09.5V 


60 [30-110] 


lO^-lO^ 


E-O 


(4) 


rho Oph 


136!| 


22,000 K, B2V 


250 [180-340] 


1-5 lO** 


E-0,C 


(5) 


N2023N 


472+J§ 


23,000 K,B 1.5V 


550 [100-1370] 


lo-'-io^ 


E-0,C 


(6) 



° Estim ated by Hipparcos iPerrvman et al.L[T997h and lHooeerwerf et alj l l2000h . * Incident FUV radiation field expressed in units x of the lDrain^ 
( fl978[) average interstellar radiation field with uncertainties in between brackets ilXmin, Xmaj])- Proton gas density n = n^jo +2 riu^ derived from 
different observational constraints. Geometry codes : S-C : Spherical Cloud; E-O : Edge-On; C : Cavity. References : (1) Habart et al. (20011); 
(2) Boulanger et al. ( 1988); (3) Rapacioh et al. (2006); Berne et al. (2007); (4) AbergeUt al. (2()03); Tevssie r et aL (2004) ; .Habart et aL C2005i) ; 
|Petvet aljr2005ll ; (5) lHabart et alJ d2003h ; lMotte et al.l ( ll998h ; r6) lField et alJ ( ll998l) ; fcomDi¥gne et alj ( l2008ll 

Table 2. H2 observed line intensities, populations in each upper level, rotational excitation temperatures, ortho-to-para ratio, and 
column densities. 



Source 


■"O-O i'(0) 

(erg s"' cm"^ sr"') 


/" 

0-0 

(erg s"' cm"^ sr"') 


■"O-O S(2) 

(erg s"' cm"^ sr"') 


■"O-O S{3) 

(erg s"' cm"^ sr"') 


■"l-O S(l) 

(erg s"' cm"^ sr"') 


Ref. 


L1721 

California 

N7023E 

Horsehead 

rho Oph 

N2023N 


6.9e-06[1.5e-06] 
4.6e-06[7.0e-07] 
3.0e-05[1.3e-06] 
1.2e-05[2.8e-06] 
6.6e-05[1.4e-05] 
5.5e-05[1.2e-05] 


2.3e-05[l.le-06] 
1.3e-05[1.3e-06] 
1.5e-04[l.le-05] 
4.3e-05[9.4e-06] 
2.5e-04[5.5e-05] 
1.5e-04[3.2e-05] 


9.1e-06[2.6e-06] 
6.4e-06[3.2e-06] 
1.5e-04[2.0e-05] 
5.4e-05[1.3e-05] 
2.4e-04[4.8e-05] 
l.le-04[4.0e-05] 


<1.0e-05 

5.5e-06[2.7e-06] 

9.2e-05[4.6e-05] 

5.9e-05[1.3e-05] 

8.6e-05[8.6e-06] 

6.2e-05[3.1e-05] 


7.5e-06[1.5e-06] 
2.3e-05[5.7e-06] 
l.le-05[2.2e-06] 


(1) 

(1) 

(1) 

(1,2) 

(1,2) 

(1,2) 



Source 


(0,2) 

(cm-2) 


(0,3) 

(cm-2) 


A?* 
(0,4) 

(cm-2) 


\0.5) 

{cm.--) 


"(1,3) 

(cm-2) 


^(0,y>4) 

^(0,J<4) 


L1721 

California 

N7023E 

Horsehead 

rho Oph 

N2023N 


4.2e+19[9.4e+18] 
2.8e-l-19[4.2e+18] 
1.8e+20[8.1e+18] 
7.1e+19[1.7e+19] 
4.0e+20[8.6e+19] 
3.3e+20[7.3e+19] 


5.2e+18[2.5e+17] 
3.0e+18[2.9e+17] 
3.4e+19[2.6e+18] 
9.8e+18[2.1e+18] 
5.8e+19[1.2e+19] 
3.5e+19[7.2e+18] 


2.6e+17[7.3e-l-16] 
1.8e+17[9.1e-Hl6] 
4.2e-Fl8[5.7e-Hl7] 
1.5e+18[3.8e-l-17] 
6.9e+18[1.3e+18] 
3.0e+18[l.le+18] 


<6.2e+16 

3.4e+16[1.7e+16] 

5.7e+17[2.8e+17] 

3.7e+17[8.3e+16] 

5.4e+17[5.4e+16] 

3.8e+17[1.9e+17] 


1.0e+15[2.0e+14] 
3.0e+15[7.7e+14] 
1.9e+15[3.7e+14] 


0.007[0.002] 
0.007[0.004] 
0.022[0.003] 
0.023 [0.008] 
0.016[0.005] 
0.009[0.004] 



Source 


^ 42 


53 


^ 64 




OTP^, 


OTP^, 




(K) 


(K) 


(K) 


(K) 






L1721 


206.[ 19.] 


<305. 










California 


208.[ 25.] 


301. [38.] 






2.3 




N7023E 


269. [ 11.] 


327.[ 43.] 






1.3 




Horsehead 


264. [ 30.] 


398.[ 49.] 


458.[ 92.] 


747. [208.] 


1.5 


2.7 


rho Oph 


252.[ 22.] 


290.[ 18.] 


327.[ 20.] 


543.[48.] 


1.1 


2.6 


N2023N 


221. [26.] 


300.[ 45.] 






1.6 





Source 


N(H2) 

(cm--) 


N(H2) I, 
(cm-2) 


N(H2^)64 

(cm -) 


N(H2) ^5 
(cm-2) 


L1721 


4.9e+20 


5.1e+19 






California 


3.2e+20 


3.0e+19 






N7023E 


1.5e-t-21 


2.9e-H20 






Horsehead 


6.1e-^20 


5.7e-Hl9 


7.5e+19 


6.0e+18 


rho Oph 


3.6e+21 


6.4e+20 


l.le-^21 


2.3e+19 


N2023N 


3.5e+21 


3.5e+20 







" Intensities at the emission peak with uncertainty (Icr) between brackets. * The level populations A'(y„_y) are derived under the assumption that 
the line is optically thin. Excitation rotational temperature, Ty^j, , derived from the ratio of the A'(o,./2) and A'(o,y|) column densities. Fit of the 
ortho-to-para (OTP) ratio of the Ji - J2 levels and assuming thermalized emission at the temperature Ty^j, . ' The total H2 column density inferred 
from the A'(o,j,) column den sities assuming optic ally thin, thermalized emission at the temperature Ty^j, . Ref . (1) Spitzer Space Telescope; (2) ESQ 
3.6-m telescope at La Silla dHabart et al.Ll2003l) and New Technology Telescope with the SOFI instrument dHabart et am200^ ; ICompiegne et all 
1200^) . 
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Fig. 1. Average IRS spectra obtained with the low spectral resolution mode. The average areas used are marked in Fig. \3] For 
LI 721 and California, all pixels are averaged in order to increase the signal-to-noise ratio, since these extended PDRs are fainter 
and do not present any detected spatial variation within the field observed with IRS. In the case of California, N7023E and rho 
Oph, only the LL IRS submodule data are available. These three objects have already been observed with the CVF of ISOCAM. 
We present here the combined ISO-CVF ( dashed lines) and IRS-LL ( solid lines) spectra. The IRS wavelength coverage allows us to 
detect several strong H2 pure rotational lines from 0-0 S(0)-S(3) at 28.2, 17.03, 12.29, and 9.66 jim, the aromatic band features at 
6.2, 7.7, 8.6, 11.3 jxm, the dust mid-IR continuum emission and the fine structure lines of ionized gas [Nell] at 12.8 jxm, [SIII] at 
18.7 and 33.4 pm and [Sill] at 34.9 jim. The dashed vertical lines delineate the H2 0-0 S(0)-S(3) lines position. The dotted vertical 
lines delineate the aromatic band feature position. 



maps were obtained using the "spectral mapping" mode with a 
number and size of steps suitable for each PDRs in order to cover 
the infrared emitting region. Projection onto the same spatial 
grid (2.5"/pixel) and convolution into the same beam (Gaussian 
with a full width high maximum of 10.7") were performed to 
compare the different line maps obtained. In the high spec- 



tral resolution mode, several individual sky positions or smaller 
spectral line maps were obtained. We used the S 17 version of the 
data delivered by the Spitzer Science Cente r. For details of our 
obs ervations and data reduc tion method see iJoblin etal] (l2005h 
and lCompiegne et al.1 (l2007h . 
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Fig. 2. Average IRS spectra obtained with the high spectral resolution mode for three objects of our sample ( LI 721, California, and 
N7023E). In each panel, the line spectrum is shown with the local rms noise (for clarity, an arbitrary offset value has been added to 
the noise). The black and green lines represent the line spectrum and the noise, respectively. With increasing wavelength, the panels 
show H2 pure rotational lines at 28.2, 17.03, and 12.29 /jm. 



We present in Fig.[T]the averaged IRS spectra obtained with 
the low spectral resolution mode in the short-low (SL, 5.2 - 14.5 
jum) and long-low (LL, 14 - 38 fim) modules. The spectra are 
averaged in a box marked in Fig. [3] showing maps in the H2 0- 
S(l) line emission at 17.03 fim. The averaged area centered at 
the emission peak covers the emission zone observed by both the 
IRS SL and LL submodules. For L1721 and California, all pix- 
els are averaged to increase the signal-to-noise ratio, since these 
extended PDRs are fainter and do not present any detected spa- 
tial variation within the field observed with IRS. In the case of 
California, N7023E and rho Oph, only the LL IRS submodule is 



available. These three objects have already been observed with 
the Circular Variable Filter (CVF A/AA -45, 5.15-16.5 yum) of 
ISOCAM. We present the combined ISO-CVF and IRS-LL spec- 
tra in Fig.[T] 

The IRS wavelength coverage allowed us to detect several 
strong H2 pure rotational lines from 0-0 S(0) to S(3) at 28.2, 
17.03, 12.29, and 9.66 /im, the aromatic band features at 6.2, 7.7, 
8.6, 11.3 ;um, the dust mid-IR continuum emission and the fine 
structure lines of ionized gas [Nell] at 12.8 /im, [SIII] at 18.7 and 
33.4 fim and [Sill] at 34.9 fim. The rotational lines 0-0 S(4), S(5) 
at 8.03 and 6.91 fim, blended with the strong aromatic feature at 
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N7023E 

H2 0-0 S(1) 



10 erg s cm sr 
. 13.0 



Iff ff'l 

3ff'\ 
S ff'\ 
30" I 

S ff'\ 



1 

I 



J2000R.A. 



rho_Oph 

H2 0-0 S(1) 



1 0"^ erg s"^ cm"'* sr"^ 
_ 31.0 




-24° 2ff 3ff 



76* 25" 58^56" 54' 52" 5(f 4^ 
J2000 R. A. 



Horsehead 

0-0 S(1) 



10'^ erg s ' cm'^ sr ' 



ff'l 



27 30" 



28 3ff' 



-2° 23 ff' 



J 



^4ff°5^Sf 5ff 55" 54' 53^ 
J2mR.A. 



N2023N 



H2 0-0 S(1) 



1 erg s cm sr 
. 21.0 



H'30" 




-2° ^3 30' 



^4r4ff 3f 36= 3€ 3t 
J2000 R. A. 



Fig. 3. Map of the H2 0-0 S(l) emission line obtained with the IRS low spectral resolution mode. The black boxes show the region 
where the spectra of Fig. \l\have been averaged. Figures of map in the other H2 0-0 S( 0), S(2 j, and S(3) emission lines are shown in 
the Appendix. 



1.1 yum and possibly the ionized argon line [Aril] at 6.98 ^m, are 
marginally detected in some sources. In the CVF spectra of the 
CaUfornia, N7023E and rho Oph, the rotational line 0-0 S(2) at 
12.29 pm appears blended with the [Nell] line at 12.8 pva be- 
cause of the CVF lower spectral resolution. To estimate the flux 
in the 0-0 S(2) line for these objects, we used the high spectral 
resolution IRS data and/or previous ISO observations (see be- 
low). In the case of the faint L1721 PDR, the 0-0 S(2) and S(3) 
lines are marginally detected in the IRS low spectral resolution 
data. We used the high spectral resolution data to estimate the 
flux in the 0-0 S(2) Hne (as described below). For the 0-0 S(3) 
line, an upper limit was derived from the low spectral resolution 
data. 

For L1721, California and N7023E, we present in Fig.|2]the 
spectra obtained with the high spectral resolution mode in the 
short-high (SH, 9.9-19.6 yum) and long-high (LH, 18.7-37.2 pm) 
modules. The data are averaged in the module slits. The wave- 



length coverage allows us to detect the 0-0 S(0), S(l) and S(2) 
Unes at 28.2, 17.03, and 12.29 pm. Both the S(l) and S(0) lines, 
which are visible in the low spectral resolution data, are clearly 
detected in the high spectral resolution data. We find good agree- 
ment between the low and high spectral resolution measure- 
ments. The S(2) line is clearly better detected in the high spec- 
tral resolution data. Nevertheless, for one object (California), the 
S(2) line is marginally detected with only 2cr detection in ampli- 
tude. 

When available, we complement and compare our mea- 
surements with ISO short wavelength spectrometer (SWS, 
/i/A/l ~ 1500-2500) data. For one object of our sample (rho 
Oph), a series of pure rotational lines (0-0 S(0) through S(5)) 
have been detected using SWS at several positions towards the 
infrared emitting region (Habart et al., 2003). The SW S posi- 
tions a re localized in the IRS-LL field (see Fig. 1 in lHabart et al.l 
(l2003h and Fig.[3]in this paper). For the first rotational S(0) and 
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S(l) lines detected with both IRS/Spitzer and SWS/ISO, convo- 
lution into the same beam (SWS beam -20") were performed 
in order to compare the data. We find good agreement between 
the IRS and SWS data with differences on the lines flux <30%. 
To estimate the flux in the other rotational lines S(2)-S(5), we 
spatially averaged the SWS positions data situated over the box 
marked in Fig. [3] 

In Table|2l we report the observed intensities of the first four 
H2 pure rotational lines. Line fluxes and error bars were gener- 
ally obtained by integrating the average spectra within the lines. 
Typical error bars are about 10-30%. In Table|2] we also report 
the intensities for the H2 1-0 S(l) rovibrational line obtained 
from high spatial resolution ground based observations when 
available. The 1-0 S(l) data have been spatially averaged over 
the box marked in Fig. [3] 

3.2. Column densities, temperatures, and ortlio-to-para ratio 

The intensity /„/ of a transition between two levels u (upper) and 
/ (lower) can be converted into the population A^„ of the upper 
level, under the assumption that the line is optically thin so that 
lui = NuAuihv/4jT. The level populations A^„ derived for each 
transition are given in Table |2l We did not correct for extinc- 
tion for the mid-IR H2 rotational line intensities, as the effect 
of reddening is within our error bars for Av=5-15, which is the 
typical range of values for our PDRs (see references in Table[TJ. 
It is only at shorter wavelengths that dust extinction will signif- 
icantly attenuate the emission. For the 1-0 S(l) line at 2.12 /im, 
we correct for dust attenuation using a factor A/f /Ay - 0.14. 

Based upon the Nu column densities, we present the excita- 
tion diagrams in Fig.|4] In Table|2] we report the rotational exci- 
tation temperatures, Ty^y,, derived from the ratio of the A^(0, J2) 
and N{0,Ji) column densities. Because the ortho-to-para ratio 
could show departures from equilibrium, the quantity Ty^y, is es- 
timated when J2 and J[ are either both odd or both even. In Table 
|2] we give the rotational excitation temperatures T42, ^53, Tm, 
and Tt5. We find that the rotational diagrams exhibit a noticeable 
curvature. A single temperature cannot describe the full set of 
observed line intensities and the inequalities 7^42 < Ts3 < Tm < 
Tys are a manifestation of the curvature apparent in the rotational 
diagrams. The v=0 rotational distributions are characterized by 
an excitation temperature between ~200 and 700 K. In order to 
fit the rotational diagrams, a combination of at least two H2 gas 
components is required: one cool/warm (~100-300K) and an- 
other warm/hot (~300-700K) with much lower column densities 
(about a few percent of the main cool/warm component). 

In Table |2] we also report a fit of the ortho-to-para ratio 
(OTP) for the Ji to J2 levels with 7i = 3, ^2 = 5 and 7, = 5, 
J2 - 1 and assuming thermalized emission at the temperature 
T j^j^ . Unfortunately, not enough H2 lines have been measured 
for LI 721 to do a proper fit of the OTP ratio. For California and 
N7023E, the OTP ratios derived are about -1-2 for the /i = 3 
to ^2 = 5 levels. For the rotational levels Ji - 5 io J2 - 7, not 
enough H2 lines have been measured in the CVF spectra to fit the 
OTP ratio. In the moderate excited PDRs, such as Horsehead and 
rho Oph, the OTP ratios derived are about ~ 1 for the J\ = 3 to 
^2 = 5 levels, while about ~3 for the J\ = 5 to ^2 = 7 levels. The 
local thermodynamic equilibrium ratio OTPlte is equal to 3 at 
the relevant temperatures {T >200K). The non-equilibrium be- 
havior has already been noted in the ISO ob servations of rho Oph 
(Habart etal.', '2003) and other PDRs (e.g.. lMoutou et alill999t 
[puente et al., 1999). The H2 column densities inferred from the 
measured line intensities and assuming thermalized emission at 
the temperatures derived are also reported in Table |2] 



In spite of their simple appearance, any interpretation of 
these diagrams has to take several factors into account. First, 
the gas temperature could vary very rapidly through the PDR 
layer, and it is likely that with our spatial resolution (~10" or 
~0.02 pc at 400 pc) the temperature structure is not spatially re- 
solved. Second, only the lowest rotational levels are collisionally 
excited, and this holds up to different J for low- and high-FUV 
field and gas density PDRs. Third, the excited rotational popu- 
lations can be affected by UV pumping. An analysis using PDR 
models is thus necessary. 

4. PDR model 

4.1. Description 

To analyze the H2 emission observations, we us ed an updated 
version of the Meudon PDR code described in iLe Petit et al.l 
(12006) . The model considers a stationary plane-parallel slab of 
gas and dust illuminated by an ultraviolet radiation field coming 
from one or both sides of the cloud. It solves the chemical and 
thermal balances iteratively at each point of the cloud. 

Here, we present a grid of isochoric models with radiation 
field coming from one side. The adopted densities are riu = 10^, 
5 10^, 10"^, 5 10^, and 10^ cm"^ in the range of values prevail- 
ing in our PDR sample (see Table [TJ. For the incident radiation 
field, we use the Draine UV radiation field scaled by the factor 
X. The adopted scaling factors of the Draine radiation field at the 
cloud edge sxt x - 0.5, 5, 50, 500, 5000 in the range of values 
prevailing in our PDR sample (see Table[T]i. 

For the gas-phase abundance, our assumed abundances are 
consistent with estimates for the average abundances in the in- 
terstellar diffuse and tra nslucent medium as derived from UV ab- 
sorpti on me asurements (ISavage & SembachLll996tlMever et all 
199"^ [T998I) . For the grain properties, the size distribution is 
assumed to be a power-law function following the MRN law; 
the UV extinction curve is taken as equal to a standard galac- 
tic expression ; the gr a in absorption cross sec tions is taken from 
iDraine & Le"3 (|1984|) ; iLaor & DrainS (Il993l) . These values ai-e 
determined for spherical particles of graphite and silicates for a 
sample of radii greater than 1 nm. Values corresponding to radii 
smaller than 50 nm are considered to mimic polycyclic aromatic 
hydrocarbon (PAH) properties. The heating rate due to the pho- 
toelectric effect on small dust grains (including PAHs) is derived 
from the formalism of Bakes & Tielens (1994) with some sig- 
nificant upgrades (see Le P etit et all 120061) . A summary of the 
most important model input parameters are given in Table[3] 

4.2. H2 specifications 

For the coUisional excitation and de-excitation of H2, we 
adopt h ere the quantu m mechanical calculations given by 
Le Bourl ot et al ] (Il999h . It must be underlined that the cross 
sections for collisional excitation and deexcitation of H2 are 
still not completly known. The quantum mechanical calculations 
adopted here contain only rate coefficients for the first fifty levels 
(T <20,000 K) and do not include reactive transitions. However, 
the populations of the (v=0, J < 10) levels, where collision data 
exist, shoul d be only secondarily affected by the missing col- 
Hsion data dShaw et al.L 120051; iKaufman et all [2006 ). Only for 
high excited levels (e.g., the 1-0 S(l) line) and high densities 
and UV fields (e.g., Orion Bar for which the population of the 
excited states is determined either by collisions or by collision- 
ally modified FUV pumping) will the rate uncertainties affect the 
emission line intensities. 
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Table 3. PDR model parameters 



Parameter 


Value 


Comment 




5 10-" 


cosmic ray ionization 


Vturb 


2 km 


turbulent velocity 


Extinction 


standard galactic 


parameterization of Fitzoatrick & Massa f 1990) 


Rv 


3.10 


Av / E(B-V) 


a) 


4.20e-01 


dust albedo 


g 


6.00e-01 


dust anisotropy factor 


G 


l.OOe-02 


Dust to gas mass ratio 


P 


3 g/cm3 


Dust density (g/cm3) 


a 


3.50e+00 


MRN dust size distribution index 




l.OOe-07 cm 


dust minimum radius 




1.50e-05 cm 


dust maximum radius 


Rf 


1.5 10-'« cm' S-' 


H9 formation rate (Habart et al., 2004) 


qH-H2 




H + H9 collision rate fLe Bourlot et al.. 1999) 


He/H 


0.1 


O/H 


3.2 10-" 


Gas abundance (Mever et al.. 1998) 


C/H 


1.3 10-" 


Gas abundance (Savage & Sembach,^996) 


N/H 


7.5 10-^ 


Gas abundance (Mever et al.. 1997) 



Although there is a consensus that H2 forms on the surface 
of dust gr ains (Gould & Salpeter, 1963; Hollenbach & Salpeter, 
ll97lHjmaL ll975: Duley & WilHams, 1984), the mechanism has 
not yet been fully understood. This is partly because we are 
ignorant of interstellar grain properties and we do not under- 
stand surface reactions in the interstellar context. Here, we 
adopt a simple empirical prescription using an H2 formation 
rate as derived from astronomical observations. In the ISM stud- 
ies, the commonly adopted value of the formation rate is 3 

10"'^ c m^ s~' inferred from UV ab sorption lines from diffuse 
clouds (IjuraL 1 1 975t iGrv et all I2OO2I) . However, studies based on 
ISO observations suggest that a higher rate of about five times 
this standard value is required to account for H2 emission in 
moderately excited PDRs (such as rho Oph, S140, and IC 63, 
[Habart et al., 2004). Such a high H2 formation rate in warm 
gas could be explained by a formation mechanism involving 
a chemisorbed H atom , analogous to the mechanism presented 
bv lCazaux & TielensI (12004). Based on this, we assume a single 
constant value of Rj - 5 x R'j. This assumption is discussed in 
Sect. ID 

Another parameter is the internal energy of the nascent H2. 
The energy released by the nascent molecule (4.5 eV) is dis- 
tributed between the grain excitation, internal energy, and ki- 
netic energy of released molecules. The branching ratios are not 
known and probably depend on conditions in the cloud and the 
nature of the grain itself. Using standard option based on an 
equipartition argument, we assume that the internal energy (1.5 
eV) is distributed with a Boltzmann law among all the energy 
levels. This assumption is discussed below. 

4.3. H2 excitation 

The excited states of H2 are populated by inelastic collisions 
with gas phase species, UV pumping or the formation process. 
The first pure rotational lines (e.g., 0-0 S(0) and S(l)) essentially 
result from collisional excitation, since their upper states (v - 0, 
J -2 and 3) are relatively low lying and their critical densities 
are low even at low temperatures («cn> < 10"* cm-"* for T > 100 
K,|Le Bourlot et al., 1999). These line intensities depend mainly 
on the gas temperature at the photodissociation front. In contrast, 
the high rotational and ro-vibrational lines (e.g., 1-0 S(l)) result 



essentially from the decay of electronically excited states that 
are pumped through the absorption of FUV photons. The line 
intensities do not reflect the gas temperature but the fraction of 
the FUV photon flux pumping H2. UV pumping could also con- 
tribute significantly to the excitation of the pure rotational 0-0 
S(2)-S(5) lines, since their upper states (v = 0, 7 =4-7) are rela- 
tively high and their critical densities are high even atmoderate 
temperatures («c„7 > 10 cm- for T < 500 K, lLe Bourlot et all 
1999) . For a given line, the relative importance of UV pump- 
ing and collisions depends on the PDR physical conditions. In 
our low n and regions (such as LI 721, California with n ~ 10^ 
cm-^,;^' ~5-20), the rotational J > 3 levels are mainly controlled 
by ultraviolet pumping rather than collisions. For moderate n and 
X regions (such as N7023E, Horsehead with n ~ 10* - 10^ cm--*, 
X ~50), collisions dominate the excitation for the 7 < 3 - 4 lev- 
els, while UV pumping contributes mostly to the excitation of 
the 7 > 4 - 5 levels. For higher n and x regions (such as rho 
Oph, N2023N with n ~ 10"* - 10^ cm-^,;^ -250-550), collisions 
dominate the excitation until the 7 = 5-6 levels. 

The internal energy of the nascent H2 can also specifically af- 
fect the level populations. However, for an equilibrium between 
photodissociation and reformation, the ratio of the formation to 
UV pumping rates is only about ~ 15/85. Of all the UV pho- 
tons absorbed by H2 only 10 to 15 % lead to dissociation. Thus, 
unless the level distribution of newly formed H2 is strongly con- 
centrated toward a small number of high energy levels, the H2 
formation excitation will not specifically affect the H2 spectrum. 



5. Comparison between PDR model and 
observations 

In this section, we compare the models predictions to the obser- 
vations. Figure |4] presents the excitation diagrams as observed 
and as predicted for each object. For each object, H2 line in- 
tensities were computed from linear interpolation of our grid of 
models with the suitable values of the FUV incident radiation 
field and gas density. Error bars due to uncertainties in the FUV 
field have been estimated. In order to compare the predicted and 
observed line intensities, the PDR is tilted with respect to the 
line of sight. For simplicity, we adopt the same viewing angle 
between the line-of-sight and the normal to the PDR for all the 
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models. We assume ^ 80°, which is suitable for PDRs viewed 
nearly edge-on as most of our PDRs are. The optically thin H2 
lines surface brightness are enhanced by a factor l/cos(0) ~6 rel- 
ative to the face-on surface brightness. 

Figure |5] shows the intensity of the rotational lines 0-0 S(0) 
through S(3), as well as for the ro-vibrational line 1-0 S(l), as 
observed and as predicted as a function of x and for different 
n. Error bars for uncertainties in the FUV field are given. The 
same viewing angle between the line-of-sight and the normal 
to the FDR is applied for all the models. Figure |5] also presents 
the kinetic temperature Tg at the H/H2 photodissociation front 
(n(H°)=n(H2)) obtained in the grid of models. At the edge of the 
cloud, the main heating gas process is the photoelectric effect 
on dust. The FUV photons heat the gas by means of photoelec- 
tric emission from grain surfaces and PAHs. As expected, Tg 
increases with;^^ since the heating photoelectric rate FpE X- 1" 
the model, both the photoelectric effect heating and UV pumping 
increase with;^', which implies strong enhance by several orders 
of magnitude of the H2 rotational and ro-vibrational emission 
line intensities (see Fig.|5]). 

Comparing the excitation diagrams as observed and as pre- 
dicted (Fig. m, we find that for all the objects - except for the 
less excited PDR LI 721 - the models account for the first low 
rotational levels 7 = 2, 3 (or the 0-0 S(0) and S(l) Hnes). These 
low levels probing the bulk of the molecular gas at moderate 
temperature (80< T <400 K) are in fact reproduced well by 
the PDR model calculations as shown in the upper panels of 
the Fig. |5] For the ro-vibrational v = 1, 7 = 3 level (or the 
1-0 S(l) line) probing the UV pumped gas, the model also fits 
the data well when available. In the Horsehead, rho Oph, and 
N2023N, the observed 1-0 S(l) line intensity is in fact in good 
agreement with the PDR model UV pumping predictions (Fig. 
ID and lower left panel in Fig. |5]|. On the other hand, the PDR 
model underestimates the column densities of the excited rota- 
tional levels with 7 = 4, 5 (or the 0-0 S(2)-S(3) Hnes) by large 
factors (factors >10) in many of our PDRs - L1721, California, 
N7023E, Horsehead - as shown in Fig. |4] and middle panels of 
Fig- El In the lowest excited L1721 PDR, the discrepancy be- 
tween model and data starts even from 7 = 3 (i.e., the 0-0 S(l) 
line, see Fig. |4]or upper right panel in Fig. |5]). From the inten- 
sity of the rotational lines 0-0 S(2) and S(3) as a function of x 
(middle panels of Fig.|5]), one can indeed clearly see that the dis- 
crepancy between models and observations mostly concern the 
low/moderate PDRs with x ^ 100. In these low/moderate UV 
environments, photon heating (by pumping and/or inelastic col- 
lisions via photoelectric effect) as described in the model cannot 
account for the observed H2 excited rotational line intensities. 
For the higher excited PDRs with > 100 (such as rho Oph, 
N2023N), we find conversely that the model roughly reproduces 
most of the observed excited rotational H2 lines. However, for 
N2023N, it is difficult to conclude about the uncertainty on;^'. 

The discrepancy between the model and the observations 
cannot be due to uncertainties on the PDR inclination, since the 
limb brightening factor l/cos(0) can be at most a few only of 
what we assumed. Furthermore, enhancement geometrical fac- 
tors would increase all the H2 lines intensities, producing too 
much emission of the low -J rotational and rovibrational lines and 
not fitting the H2 line intensity ratios. Considering the error bars 
for uncertainties in the FUV field shown in Figs.|4]and|5] we also 
find that the discrepancy between the model and the data cannot 
come from uncertainties on the strength of the FUV incident ra- 
diation field X- Similarly, uncertainties on the local gas density 
n cannot be invoked as the predicted intensities for excited ro- 
tational lines are comparable for models with n = 10^ - 10^' 



cm"^ (see Fig.|5]for;i' < 100). Density structure changes would 
also have minor effects, since modeling the H2 excitation mainly 
depends upon the average gas density in the H2 emission zone. 
Detailed calculations that take a steep gas density gradient into 
account as observed fro m the Horsehead PDR edge from previ- 
ous complementary data jHabart et al.L l2005l) led in fact to sim- 
ilar results. However, dumpiness that can lead to a substantial 
reduction of the effective optical depths could increase the por- 
tion of hot gas and explain a significant amount of rotationally 
excited H2. 

6. Comparison with H2 observations of the diffuse 
interstellar medium in the Solar Neighborhood 

Early far-UV absorption studies with the Copernicus satel- 
lite, as well as far-UV spectroscopic data obtained more 
recently with FUSE show significant amoun t s of excited 
H2 in diffuse clou ds (Spitzer & Cochran', '1973'; 'Spitze r et al. , 
'19741; ISpitzer & Morton,' 1976; Savage et al., 1977; Gr v et al. . 
2002; Browning et al., 2003; Gillmon et al., 2006). In particular, 
,Grv et al.i (i2002) analyze FUSE H2 spectra over three lines of 
sight in the Chamaeleon clouds probing gas heated by the mean 
Solar Neighborhood radiation field. These data, which include 
the populations of the lowest two 7=0 and 7=1 levels, clearly 
reveal a multi-temperature structure. The 7 > 3 levels charac- 
terized by an excitation temperature of T52 ~200-250 K lie far 
above the temperature derived from the lowest level transitions 
with TiQ -60 K. Using the M eudon PDR code withx - 0.4 and 
n ~ 80 cm"^, iNehme et al.i t2008) find a good fit to a number 
of observables, including the column density of H2 in the 7=0 
and 1 levels, but their model fails to account for the H2 in levels 
7 > 3. The discrepancy between the observed and predicted col- 
umn densities of these levels is about one order of magnitude, a 
difference similar to that we report in low-excitation PDRs. 

The H2 excitation in PDRs and in the diffuse interstellar 
medium differ in the power radiated per hydrogen atom. The 
power radiated from the first rotational levels per nucleon, i.e. 
Pff. ~ Iff./Nn X An (in erg s"' H"'), where Ih' is the inten- 
sity of the high rotational lines and Nh the total hydrogen col- 
umn density of the emission zone, is much stronger in low- 
excitation PDRs than in diffuse clouds. For the L1721 PDR with 
/h- ~ /()-{) + h)-o s(2) ~ 4.8 10"^ erg s"' cm"^ sr"' and A^^ ~ 
2xNh2 ~ 10^' cm-2 (see Table |2]i, we find Ph-^ ~ 6 lO^^s gj.g 
s ' H ' . For the diffuse interstellar medium towards HD 102065, 
the observed H2 column densities translate into Ih^ ~10"^ erg 
s ' cm"^ sr adding emission in the S(l) and S(2) lines. With 
Nh ~ lO^i cm-2, we find Ph- ~ 2 lO^^s erg s"' H \ This 
value is more than one order of magnitude lower than the LI 721 
value. The difference is even greater for other PDRs in our sam- 
ple. This result indicates that the energy deposited per hydro- 
gen atom is much higher in PDRs than in diffuse interstellar 
clouds. This is not surprising since physical conditions in PDRs 
differ from those in the diffuse interstellar medium. In particu- 
lar, the gas kinetic energy per nucleon and its dissipation rate 
could be higher, because PDRs in our sample are all in massive 
star-forming regions where turbulence may be powered by the 
mechanical action of luminous stars on their environment. Thus, 
the mechanical processes propose d to account for the H2 excita- 
tion in diffuse interstellar cl ouds jFlower & Pineau Des Foreti, 
119981: iFalgaroiie et al.L I2b05h could also contribute to the H2 ex- 
citation in PDRs, even if the power radiated per H atom in the 
rotational excited lines is much greater. This possibility is further 
discussed in Sect. 17.41 
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Fig. 4. Excitation diagram of Hi- Nu is the column density of the transition upper level, gu is the degeneracy of the upper level and 
T,, is the upper level energy in Kelvin. The squares show the data. The H2 ortho-to-para ratio is taken equal to 3. The error bars 
are small in size compared to the symbol size. The arrows indicate upper limits. The stars show the PDR model predictions. For 
each object, H2 line intensities were computed from linear interpolation of our grid of models with the suitable values of the FUV 
incident radiation field and gas density (see Table\I}. We use x— 4.5 and n — cm^^ for L1721; x—^^ '^^d n—5 XQ^ cm^^ for the 
California; X— 5 1 and n — 10'^ cm^^ for N7023E; x—60 andn — 10'^ cm^^ for the Horsehead; x-250 and n = IQ'^ cm~^ for rho Oph; 
X=550 and n = 5 10^ cm^^ for N2023N. Error bars due to uncertainties in the UV field are given. For comparison of the predicted 
and observed line intensities, in the model the PDRs have been tilted with respect to the line of sight with cos(0)—O.15. 
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X (Draine units) X (Draine units) 

Fig.5. Intensites of several pure rotational lines 0-0 S(0) through S(3) as well as for the ro-vibrational line 1-0 S(l) predicted by 
the PDR model for different gas densities - n — 10^ cm^^ ( dashed lines), n = 10"* cm"-' ( solid lines), n — arT^ ( dotted lines) - as 
a function of the strength of the FUV incident radiation field, X- The kinetic temperature at the photodissociation front Tg obtained 
in the grid of models is also presented in the lower right panel. The blue points with error bars show the observed intensity for our 
PDR sample. From left to right, we have L1721 (x-4.5), California (x-19), N7023E (x-51 ), Horsehead (x-60), rho Oph (x-250), 
and N2023N (x—550). Horizontal error bars due to uncertainties in the FUV field are given. Vertical error bars for uncertainties in 
the intensities are also given. The arrows indicate the upper limits. For comparison of the predicted and observed line intensities, 
in the model the PDRs have been tilted with respect to the line of sight with cos(6)—0.15. 



7. Origin of the rotationally excited H2 

Our Spitzer results show the presence of significant rotationally 
excited H2 in the low/moderate excited PDRs that is not pre- 
dicted by PDR models. The cause of this could be that, in the 



models, the column density of H2 is too low in the outer zones 
where H2 is excited or alternatively that a small fraction of the 
H2 gas is hotter than predicted by models. The effect should be 
selective by only significantly increasing the intensities of the 
excited H2 rotational lines (e.g., 0-0 S(2)-S(5)). The intensities 
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of the low rotational lines (e.g., 0-0 S(0)-S(1)) probing the bulk 
of the molecular H2 gas at moderate temperature, as well as of 
the rovibrational line (e.g., 1-0 S(l)) probing the UV pumping 
excitation, should mostly not be affected in order to fit the data. 
In the following, we examine some of the issues with this in 
mind. 

7.1. H2 formation rate 

The H2 line intensities depend mainly on the H2 formation rate 
Rf since it controls the location of the H/H2 transition zone in 
PDRs. The line intensities change via the variation in both the 
gas temperature Tg and the flux pumping H2 in the H/H2 transi- 
tion zone. Enhanced Rf could result in a strong increase in the 
intensity of the excited rotational lines. Previous ISO studies of 
moderately excited PDRs (such as rho Oph, S140, and IC 63) 
suggest that a high R j about five times the standard value is re - 
quired to reproduce the H2 rotational lines (iHabart et al.Ll2004l) . 
However, we find here that, even if the H2 formation rate is sig- 
nificantly enhanced, the model intensities fall s hort by an order 
of magnitude. The same conclusion is r cached bv lLi et all (|2002|) 
using ISO data of the low-excited PDR the S140 periphery re- 
gions (with ;r -10-30 and n ~ 10^ cm^^). This can be understood 
as follows. In this low/moderately excited PDR regime with low 
Xin ratios ixln < 0.01), H2 self-shields efficiently enough so 
that the H°/H2 transition zone is closer to the edge (A^ < 10^' 
cm"^), even without an enhanced H2 formation rate. H2 already 
competes effectively with dust in absorbing FUV photons and 
the gas temperature in the H'VH2 transition zone is comparable 
to the PDR edge temperature, so although the H2 line intensities 
are increased in the high H2 formation case, it does not increase 
much. We conclude that a global enhancement in the H2 forma- 
tion rates improves the comparison with observations, but is not 
sufficient to match the excited rotational lines as observed in our 
low-excited PDRs. 

7.2. Advected H2 

Our model results are based on the assumption that the PDR 
is in equilibrium. In reality, advection of molecular gas from 
the shielded cloud interior to the warm surface - created by 
turbulent motions and/or fast progression of the photodissoci- 
ation front into the cloud - can create relatively high nonequi- 
librium abundances of H2 in the outer layers. This could pro- 
duce extr a column densiti es of excited H2. However, as dis- 
cussed in IStorzer "&Hollenbachl (ll998h . there is a range of x/n 
(0.01 < x/f^ 0.2v with V the gas velocity in km/s) for which 
the effects of advection will have the strongest influence on line 
intensities compared with equilibrium values. For high;\f/n, the 
dissociation timescales are so short that advected H2 will not 
survive in the outer layers. For very low x/n, the differences 
between equilibrium and nonequilibrium line intensities are ex- 
pected to be small, since the H/H2 front is in the equilibrium 
model near the edge at N < 10^' cm"^, where most of the ex- 
cited H2 comes from. For our low/moderately excited PDRs with 
10^ < < 10"^, the condition that advected H2 will survive 
in the outer layers without being dissociated is verified, even for 
low velocity. However, the line intensities are expected not to in- 
crease significantly, since the H/H2 front is near the edge, even 
without advection. Out-of-equilibrium PDR model calculations 
by Storzer & Hollenbach (1998) predict that the integrated H2 
line intensities do not vary by more than factors of 3. However, 
their calculations concern the high-density, high-;^' case, which 



does not apply here. Moreover, local increase in the H2 lines 
emissivity at the PDR edge can be important. Further detailed 
calculations for the low-;^' case would be needed to fully test the 
non-equilibrium effects. 

Another point that one can examine is the OTP ratio expected 
in the advection case. If there is advection and the dynamical 
timescale {tj) is much shorter than the OTP equilibrium time 
(toTp), the OTP ratio must be out-of-equilibrium. Simple calcu- 
lations of these different timescales for the conditions prevailing 
in our PDR^ and assuming v ~1 km/s suggest that tf < tpTP 
for the low excited PDRs (such as L1721 with;if ~ 5, n ~1Q^ 
cm"^) while tf > tgrp for the outer layers of the moderate ex- 
cited PDRs (such as N7023E, Horsehead with;i' ~ 50, n -10"* 
cm"^). Consequently, the OTP ratio expected in the advection 
case should be equal to 1 or 3 according to the PDR zone stud- 
ied. Unfortunately, for many of our PDRs, not enough H2 lines 
have been measured to do a proper fit of the OTP ratio of the ex- 
cited rotational levels (see Sect. [3]). Nevertheless, for moderately 
excited PDRs (Horsehead, rho Oph), where several excited H2 
rotational lines have been measured, we derive the OTP ratio of 
about ~ 1 for the 7=3 to 5 levels while ~3 for the J -5 to 7 lev- 
els. The different OTP findings could result from the fact these 
rotational levels probe different PDR zones where the OTP con- 
version is more or less efficient. In either case, we suggest that 
these findings should not eliminate advection as a solution to 
extra warm/hot H2. However, more sensitive searches to better 
constrain the OTP ratio of the excited H2 gas, as well as con- 
straints on the kinematics in the H2 emission zone, are needed 
for any conclusion. 

7.3. Local heating via tiie piiotoeiectric effect 

One alternative solution is to locally increase the gas temper- 
ature in the PDR outer zone. At the edge of the cloud, the 
main heating gas process is the photon heating via photoelec- 
tric effect on dust. A local increase in the photoelectric heat- 
ing rate can have various causes. In particular, there is increas- 
ing observational evidence of evolution of the very small dust 
particles, including PAHs and very small grains (VSGs) in the 
illumina ted parts of the cl oud (Boulanger et al., 1988, 1990; 
Bernard et al., 1993; Abergel et al., 2002;" Rapacioli et aUI2005l: 
Be rne et al ., 20 07; Comp iegne et al., 2008.). Increasing the abun- 
dance of the smaller grains that dominate the p hotoelectric effect 
would lead to higher gas temperatures (e.g., iBakes & Tielensl 
1994; Weingartner & Draine, 2001). Previous ISO study of the 
low-excited PDR L1721 bv lHabart et an(l200lh shows that abun- 



' The advection time can be given as tf = No/(n x i') ~ 3 10* yrs x [ 
(llcm/s) / V ] X [ (10'* cm"^) / n] with A'o - 10"' cm"^ the total hydrogen 
column density in the surface PDR zone. The dominant mechanisms 
that control the OTP ratio conversion via gas-phase spin exchange re- 
actions are collisions with H atoms and/or protons. The rate coefficient 
for the proton reaction can be given by 3 10"'" cm'' s"' (Gerlich, 1990|). 
The rate coefficient for the reactive collisions with H atoms can be given 
by 8 10"" s^-i-xxi/T,) s"' (e.g., :Schofield, 1967). Due to the high 
activation barrier of the reaction with H atoms (~3900 K), the domi- 
nant process that controls the OTP ratio conversion is collisions with 
protons in cool gas (Tg < 100-200 K). Taking the gas temperature pro- 
files as predicted by the equilibrium model, we find that Tg <200 K for 
PDR withx ~ 5, while Tg <350 K for PDR with x ~ 50. Thus, for 
low-excited PDRs the proton reaction dominate all the PDR, while for 
moderatly excited PDRs, the collisions with H atoms are expected to be 
important at the PDR edge. Within the range of conditions prevailing in 
our PDRs, we estimate torp ~ 10*- 10' yrs for the proton reaction and 
toTP ~ (3 10' - 10")/;i yrs for reaction with H atoms. 
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dance variations of small dust grains should be taken into ac- 
count to reproduce the cooling gas due to fine structure lines 
([C+] 158 fim and [O"] 63 fim). Also, they predict that the H2 
lines emission is increased by a factor of ~5-10, in line with the 
present data. 

Analysis of the full IRS Spitzer data spectrum obtained in 
the frame of the SPECPDR program permits us to characterize 
the small grain populatio n. In particular, Berne et al.. (2007) and 
ICompiegne et al.l (l2008l) observed significant variations in the 
relative emission between R\Hs and VSGs in our PDRs sam- 
plfl Compiegne et al. (2008), who focus on the Horsehedad 
and N2023N, shows that excitation effects cannot account for 
these variations and interpr ets it in term of PAH/VSG relative 
abundance variations. Bern e et al. I (l2007l) extracted the mid-IR 
spectra of PAH and VSG populations. For N2023N, where spec- 
tral evolution is observed, Compiegne et al. (2008) find that the 
PAH/VSG abundance ratio is ~5 times lower in the dense, cold 
zones of the PDR than in its diffuse illuminated part where dust 
properties seem to be the same as in cirrus. For the Horsehead, 
where they did not spatially resolve any spectral variation, they 
found that this ratio is 2.4 less at the peak emission of the PDR 
than in Cirrus. However, stronger variations on spatial scales be- 
low the available Spitzer resolution cannot be excluded. To quan- 
tify this impact, the absolute grain abundance profile of each dust 
population as a function of depth should be implemented in the 
PDR model. Towards our PDRs where spectral evolution is ob- 
served, more detailed studies of the spatial emission profiles of 
the dust components, in relation with the H2 lines, are currently 
in progress (.Joblin et aL.201 1.) . 

7.4. Heating by dissipation of meclianical energy 

Here, we consider H2 heating associated with the dissipation of 
the gas kinetic energy. PDRs in our sample are all located in mas- 
sive star-forming regions, where the expansion of HIl gas, stel- 
lar winds and/or supernovae explosions are significant sources 
of mechanical energy. The energy released in shocks or by the 
dissipation of turbulence could increase the H2(7) populations 
in comparison to PDR model predictions. We find that the MHD 
shock model of Flower & Pineau Des Forets ( 19980 fits the ex- 
citation well as do as the absolute intensities of the H2 lines ob- 
served in the low-excited PDR LI 721, for a gas density n ~ lO^* 
cm"-' and a shock velocity v ~ 15 km/s. 

The dissipation rate of turbulence (Pdiss) is often related 
to the amplitude of turbulent motions (the gas velocity dis- 
persion Vrms) and the size (/,) at which the energy is injected. 
Quantitatively we have Pdiss/Mn ~ 3/2 x v^^Jli = 6 x 10"^^ x 
(vVm.s/20 km s^^f" IHijlQ pc) erg s"' H"'. This estimate quanti- 
fies the magnitude of the turbulent gas velocity, which is neces- 
sary to account for the H2 S(l) and S(2) line emission in L1721. 
High spectral resolution observations of H2 (R - Aj t^A greater 
than 20,000 achievable with the CRJRES spectrograph on the 
VLT) are needed in order to test this possibility. High spec- 
tral resolution observations of the C^ far-IR line emission from 
PDRs, to be performed with the Herschel Space Telescope, will 



N7023E and rho Oph spectra were analyzed bv [ Beme et aP l l2007h : 
Horsehead and N2023N spectra bv'Comp iegne etal.l ( 120081) . L1721 and 
California do not present any detected spatial variation within the small 
field observed with IRS. 

^ The MHD shock models have been updated, i n particular with 
the ne w cross sections of collisional excitation of H2 iLe Bourlot et alj 
jl99gh . 



also quantify turbulence in the outer PDR layers where CO is 
photodissociated. 

8. A diffuse contribution to tlie H2 emission from 
galaxies 

Spitzer spectroscopic observations have provided detections 
of the H2 rotational lines for many galaxies, in parti cular a 
large fraction of galaxies in the SINGS sample (Rousse fet all 
2007). The H2 emission of star-forming galaxies is tightly cor- 
related with the PAH emission. It is particularly striking that the 
H2/PAH emission ratio is observed to be independent of the in- 
tensity of the dust heating, as measured by the 24 yum to far- 
IR dust emission ratio. Roussel et al, (2007) argue that the H2 
and dust emission are powered by a common energy source, the 
UV light from young stars. The detection of H2 line emission 
in low-excitation PDRs, such as L1721, supports this interpreta- 
tion. Low-excitation PDRs, such as LI 721, may add up to make 
a diffuse component that may be a dominant contribution to the 
H2 and PAH emissions from star-forming galaxies. 

Rotationally excited H2 has been detected along Galactic 
lines of sight away from lunimqus star fqrmirig regions with ISO 
and Spitzer observations dFalgarone et all (20051: iHewitt et all 
2009). Thus, Galactic observations do indicate the presence of 
a diffuse component of H2 emission in galaxies. However, it is 
still unclear whether this emission comes from low-excitation 
PD Rs or if it arises from the diffuse interstellar medium. Like 
us, iFalgarone et al.l (l2005h used the Meudon code to estimate 
the contribution of low-excitation PDRs to the observed line 
emission. They concludes that the S(0) line emission may be 
accounted for by PDRs but that the model estimates fall short 
of the observed S(l) and S(2) lines by a large factor (~5-10). 
This discrepancy is similar to the one we faced in modeling our 
observations of low-excitation PDRs. The absolute line intensi- 
ties measured along their galactic line of sight are comparable 
to the intensities measured in our L1721 PDR, but the S(0)/S(1) 
line inte nsity ratio is significa ntly higher in the Galactic lines 
of sight. IFalgarone et"aL (|2005|) concludes that the UV emission 
cannot be the unique energy source and that much of the H2 line 
emission must come from the diffuse interstellar medium where 
H2 may be heated by the localized dissipation of the gas turbu- 
lent kinetic energy. This conclusion needs to be revisited in the 
light of PDR observations, but it is not necessarily contradicted 
by the correlation between H2 and PAH emissions, since the dif- 
fuse interstellar medium is also thoiight to provide a significant 
contribution to the PAH emission dbraine et alll2007i) . 

9. Conclusion 

Thanks to Spitzer spectroscopic observations, we were able to 
detect H2 pure rotational lines emission in PDRs with modest 
FUV fields and densities. The low/moderate x regions studied 
here are very widespread in galaxies and may contain a large 
fraction of the molecular gas. However, this intermediate regime 
between the diff'use cloud and the bright PDRs has been poorly 
studied. To analyze the H2 line emission observations, we used 
an updated version of the Meudon PDR code. Our results allow 
strong constraints to be placed upon the H2 excitation in the in- 
terstellar medium of galaxies. The main results from this work 
can be summarized as follows. 

The IRS wavelength coverage allows us to detect several 
strong H2 pure rotational lines from 0-0 S(0) to S(3) at 28.2, 
17.03, 12.29, and 9.66 ;um, the aromatic band features at 6.2, 
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7.7, 8.6, 1 1.3 ijm, the dust mid-IR continuum emission, and the 
fine structure lines of ionized gas [Nell] at 12.8 yum, [SIII] at 
18.7 and 33.4 fim and [Sill] at 34.9 fim. The observed mid-IR 
spectra are typical of PDRs. 

A single temperature cannot describe the full set of observed 
H2 line intensities. A combination of at least two H2 gas compo- 
nents, with one cool/warm (~100-300K) and another warm/hot 
(~300-700K) with much lower column densities (a few percent 
of the first component) is required. The ortho-to-para ratios de- 
rived are about ~1 for the first component and about ~3 for 
the second warmer component. The non-equilibrium behavior 
has already been noted in previou s ISO observations of PD Rs 
(iMoutou et al.Lll999HFuente et ali fT999: Habart etan i2003h . 

By comparing the observations with the FDR model predic- 
tions, we find that the model can account for the first low H2 
rotational line (e.g., 0-0 S(0)-S(1)) probing the bulk of the gas 
at moderate temperature, as well as the ro-vibrational line (e.g., 
1-0 S(l) observed with ground-based telescopes) probing the 
UV-pumped gas. In contrast, the model underestimates the ex- 
cited rotational fines (e.g., 0-0 S(2)-S(3)) by large factors (>10 
in some objects). In the lowest excited PDR, the discrepancy 
between the model and the data starts even from the rotational 
J - 3 level (e.g., 0-0 S(l)). This discrepancy between obser- 
vations and the PDR model predictions has the same order of 
magnitude as reported for diffuse interstellar clouds. However, 
in PDRs the power radiated per H atom in the rotational excited 
levels is more than one order of magnitude larger than in diffuse 
clouds. 

The discrepancy between the data and the PDR model could 
indicate that, in the models, the column density of H2 is too low 
in the outer zones where H2 is excited (by inelastic collisions or 
UV pumping), or, alternatively, that a small fraction of the H2 
gas is hotter than predicted by models. We discuss whether an 
enhancement in the H2 formation rate, or a local increase in pho- 
toelectric heating, as proposed for brighter PDRs in former ISO 
studies, may also work in low-excitation PDRs. An enhancement 
in the H2 formation rate improves the comparison with observa- 
tions, but the models still fall short. Further work is needed to 
quantify the impact of the evolution of very small dust particles 
across PDRs on the gas energetics. Out-of-equilibrium processes 
or mechanical heating (by weak shocks or dissipation of turbu- 
lence as proposed for diffuse interstellar clouds) are alternative 
possibilities. 

Although we cannot decide at this point, we emphasize 
the need for further development of the models. Progress is 
also expected from PDR spectroscopy of additional cooling 
lines and data on the gas kinematics, to be obtained from the 
ground and Herschel. In particular, observations of gas cool- 
ing lines of species such as at high spectral resolution with 
the Heterodyne Instrument for the Far Infrared (HIFI) will pro- 
vide missing information on the gas velocity within the PDR 
layer where CO is photodissociated. The combination of these 
data with models should help constrain the relative radiative and 
dynamical influence of stars on the physical conditions within 
PDRs. In the longer term, the Mid-InfraRed Instrument (MIRI) 
on James Webb Space Telescope (5-28 //m, diffraction-limited 
resolution of ~0.2" at 5 fim) will provide maps of the H2 rota- 
tional line emission, as well as of the small dust emission with 
unprecedented angular resolution. 

Appendix A: Map in the H2 rotational emission lines 

Figure A.l shows maps in the H2 0-0 S(0) line emission 
at 28.2 fim obtained with the IRS long-low module towards 



N7023E, the Horsehead, rho Oph, and N2023N. Maps in the 
H2 0-0 S(l) line emission at 17.03 jum are presented in Fig. [3] 
Figure A. 2 shows maps in the H2 0-0 S(2) and S(3) lines emis- 
sion at 12.29 and 9.66 fim obtained with the IRS short-low mod- 
ule towards the Horsehead. Projection onto the same spatial grid 
(2.5"/pixel) and convolution into the same beam (Gaussian with 
a full width high maximum of 10.7") were performed to com- 
pare the different line maps obtained. There is a good agreement 
in the intensity distribution between the H2 lines. The same sub- 
structures are present in the maps. The emission in the lowest 
rotational line 0-0 S(0) is somewhat more extended and peaks 
more inside the nebula. This is easily to understand since the up- 
per state of this line is relatively low-lying. The signal-to-noise 
ratio in the 0-0 S(0) line maps is lower than in the other line maps 
since this line is fainter. Towards the Horsehead, the spatial dis- 
tribution of the excited 0-0 S(2) and S(3) lines is comparable. 
The 0-0 S(3) line is sUghlty shifted nearer the edge. 
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Fig. A.l. Map of the H2 0-0 S(0) emission line obtained with the IRS low spectral resolution mode. The black boxes show the region 
where the spectra of Fig.\l\have been averaged. 
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Fig. A.2. Map of the H2 0-0 S(2) and S(3) emission line obtained with the IRS low spectral resolution mode towards the Horsehead 
nebula. The black boxes show the region where the spectra of Fig. \l\have been averaged. 
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